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Abstract. — We present a semiclassical perturbation method for the description of atomic diffraction by 
a weakly modulated potential. It proceeds in a way similar to the treatment of light diffraction by a thin 
phase grating, and consists in calculating the atomic wavefunction by means of action integrals along the 
classical trajectories of the atoms in the absence of the modulated part of the potential. The capabilities 
and the validity condition of the method are illustrated on the well-known case of atomic diffraction by a 
Gaussian standing wave. We prove that in this situation the perturbation method is equivalent to the Raman- 
Nath approximation, and we point out that the usually-considered Raman-Nath validity condition can lead 
to inaccuracies in the evaluation of the phases of the diffraction amplitudes. The method is also applied to 
the case of an evanescent wave reflection grating, and an analytical expression for the diffraction pattern at 
any incidence angle is obtained for the first time. Finally, the application of the method to other situations 
is briefly discussed. 



1. Introduction. 

The diffraction of atomic de Broglie waves by standing wave light fields [1-4] or mechanical mi- 
crostructures [5] has recently received considerable interest because of its potential application as an 
atomic beam splitter, one of the key components for the development of atom optics and interferometry 
[6]. Among the different realizations of atom diffraction gratings proposed to date, the near-resonant 
Kapitza-Dirac effect [7] (the diffraction of atom from a standing laser field) has led by far to the most 
theoretical work [8]. Widespread interest in this effect arose because the phenomenon is the quantum 
mechanical analog of diffraction of light waves by a matter grating and also because it is conceptually 
one of the simplest examples of stimulated momentum transfer between atoms and light. The theoret- 
ical description of such a transmission diffraction is greatly simplified in the Raman-Nath regime of 
diffraction [1] where the change in kinetic energy of the atoms due to diffraction is neglected compared 
to the atom-field coupling. The problem then reduces to one dimension and the diffraction grating acts 
as a thin phase grating (the interaction between the atoms and the stationary laser wave only affects 
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the phase of the atomic wavef unction). In contrast, the theoretical treatment of atomic diffraction by 
an evanescent wave reflection grating [3, 4] is known to lead to several unusual problems, which are 
related to the fact that the motion perpendicular to the diffraction grating can no longer be eliminated in 
a constant motion approximation. Indeed, the slowing down and eventually the reversal of the atomic 
motion is intrinsically associated with the reflection grating, and the problem remains necessarily two- 
dimensional. As a consequence, only a few attempts [3, 4] have been made to describe this diffraction 
process, and to our knowledge, no theoretical treatment provides a clear physical picture of atomic 
diffraction in the limit of small saturation of the atomic transition (a regime of great interest in atom 
optics experiments). 

We present in this paper a semiclassical perturbation method which permits the treatment of atomic 
diffraction by a weakly modulated potential. The method applies both for transmission and reflection 
gratings, and proceeds in a way similar to the treatment of hght diffraction by a thin phase grating. It 
is based on the evaluation of the atomic wavefunction by means of action integrals along the classical 
trajectories of the atoms calculated in the absence of the modulated part of the potential. In order to 
illustrate the capabihties and the limits of the perturbation method, we first consider the well-known 
near-resonant Kapitza-Dirac effect for which we prove the equivalence between our method and the 
Raman-Nath approximation. Furthermore, we show that the commonly-accepted validity condition 
of the latter actually leads to inaccuracies in the phases of the diffraction amphtudes. The case of 
atomic diffraction by a weakly modulated evanescent wave grating is then investigated and an analytical 
expression for the diffraction pattern at any incidence angle is derived for the first time. Finally, we 
discuss the apphcation of the method to the situation of time-modulated potentials and to the case of 
multilevel atoms. 

2. Semiclassical perturbative calculation of the diffraction spectrum. 

We present in this section the principle of the semiclassical perturbation method for calculating the 
diffraction spectrum of an atomic de Broglie wave interacting with a weakly modulated potential. For 
illustration purposes, we will here restrict the discussion to the case of a spatially modulated potential 
(the discussion of a time-modulated potential is postponed to section 5). 

2.1. Description of the model. — We consider the simple case of a two-level atom incident on 
the optical diffraction grating provided by an appropriate arrangement of laser beams. Because we 
are interested in the regime of coherent atom optics (Umit of negligible spontaneous emission), we 
restrict ourselves to the case of low saturation of the atomic transition where the reactive part of the 
atom-laser wave coupling (hght-shifts) is predominant over the dissipative part. We also assume that 
the detuning between the laser waves and the atomic frequency is properly chosen so that the atoms 
follow adiabaticaUy the optical potential associated with the light-shifted ground-state level, and that 
any Doppler effect can be neglected. The Lagrangian of the atomic system is then of the form: 

L{r,v) = Lo{r,r)-eVir) (1) 

where r stands for the position of the atomic centre of mass, and r for its velocity. In equation (1), eV{r) 
denotes the spatially modulated part of the optical potential responsible for atomic diffraction. This 
potential is assumed to be a small perturbation (e <C 1 is the perturbation parameter) compared to the 
non spatially modulated Lagrangian Lq which contains the kinetic energy term. 

2.2. Calculation of the atomic wavefunction. — 

2.2. 1 . The WKB method. — In the framework of a semiclassical (WKB) treatment of the atomic center 
of mass motion, the atomic wavefunction is evaluated by means of action integrals along classical 
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atomic trajectories. As shown in Appendix A, assuming that at time t = ti, the atoms have not yet 
interacted with the diffraction grating and that the atomic wavefunction corresponds to a plane wave of 
momentum pi, the wavefunction at position rf and time tf is given by: 



ip{r{,t[) = exp 



h 



Sir[,ti\pi,ti) 



(2) 



where: 



S{r{, t{\pu ti) = Pi • r{ti) + ' dtL[T{t), r 



v{t)\ (3) 



is the action integral along the classical trajectory v{t) solution of the Euler-Lagrange equations asso- 
ciated with the Lagrangian L (Eq.(l)), given the boundary conditions: 




(4) 

rf 

with M the atomic mass, and the requirement that r(ti) be situated outside the interaction region. Note 
that the first term in the right-hand side of equation (3) takes into account the fact that the boundary 
conditions (4) differ from the usual case where the initial and final positions of the trajectory are 
specified. This term is associated with the phase of the atomic wavefunction at the position r(ti). 

It is interesting to note that the above-described method for evaluating the atomic wavefunction 
closely resembles the way one accounts for interference or diffraction effects in conventional optics 
(where optical paths are calculated along light rays derived from Fermats principle), and therefore 
is subject to the same validity conditions. More precisely, it requires that both the amphtude of the 
wavefunction and the optical potential vary slowly on the scale of the atomic de BrogUe wavelength. It 
thus breaks down near the points where classical trajectories cross each other, e.g., near caustics or focal 
points. Another characteristics of the WKB method is that it requires the knowledge of the classical 
trajectories for the total Lagrangian L, which in practical application requires a numerical integration 
of the Euler-Lagrange equations. 



2.2.2. The perturbation method. — We show here that by taking advantage of the weakness of the spa- 
tially modulated potential £V{r), it is in fact possible to evaluate the action integral (3) perturbatively 
up to leading order in s using only the classical atomic trajectories for the unperturbed Lagrangian L 
[9, 10]. The interest of this method is that the unperturbed trajectories are often known analytically, 
and therefore allow for analytical derivations of the diffraction spectrum. 

The perturbation method proceeds as follows. We expand the action integral (3) and the actual atomic 
trajectories r{t), solutions of the Euler-Lagrange equations for the perturbed Langrangian L, in powers 
of the small parameter e: 

r(t) = roit) + sri{t) + e^r2{t) + ■ ■ ■ 
S{ri,t{\pi,ti) = So+eSi+e^S2 + --- 

Note that ro(t) corresponds to the unperturbed classical trajectories, solutions of the equations of mo- 
tion for the unperturbed Lagrangian Lq. Substituting equation (5) into (3), using equation (1) and 
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separating the different orders in e, we get a set of equations the three first of which are (see Appendix 
A): 



So = Pi-ro{ti) + dtLo[roit),roit)] (6a) 



Si = - j\tV[v^{t)\ 



(6b) 

^2 = -i^ 'dtri(t)- Vy[ro(t)] (6c) 

Equation (6a) is nothing but the exact action integral (3) in the limiting case £ = where the mod- 
ulated part of the optical potential vanishes (no atomic diffraction). As a result, even though So is 
generally known analytically, it plays no role in the characterization of the diffraction spectrum and 
therefore will not be considered in the following. More relvant is equation (6b) which describes the 
phase-shift accumulated by the atom along its unperturbed trajectory ro{t), due to the presence of the 
modulated part of the optical potential. 5*1 actually contains the first order information about the dis- 
tortion of the atomic wavefront by the diffraction grating, and can thus be used to derive the diffraction 
spectrum. This is the central point of our perturbation treatment of atomic diffraction. It is clear how- 
ever that such a method will only be accurate in the limit of small e, and we now discuss its validity 
range. Consindering equations (2) and (5), it appears that an appropriate condition is that be suf- 
ficiently small compared to h, in other words that the higher perturbation orders do not significantiy 
affect the atomic wavefunction. Using equation (6c), the validity condition of our perturbation method 
thus reads: 



dtri(i)- Vt/[ro(t)] 



ti 



< n (7) 



for all possible trajectories. In fact, it is interesting to find an upper limit for the left-hand side of 
equation (7) which permits a more transparent physical interpretation of the validity condition. Note 
first that the gradient —eW is the additional force acting on the atom due to the modulated part of 
the potential. From the classical point of view, this force is responsible for the momentum transfer 
involved in the diffraction process. As a result, the time integral of the force will be of the order of 
the maximum momentum transfer /S.pmax observed in the diffraction spectrum (see Fig. 1). Second 
|£ri(t)| corresponds to the deviation of the atom from its unperturbed trajectory. Its maximum is 
Ar^ax. the largest displacement observed after the atoms have excited the interaction region with the 
light grating (see Fig. 1). One thus obtains a condition for the validity of the perturbation method: 

^Ap„iax^r„iax < ^ (8) 

which states that the error in the atomic phase due to the integration of the classical action along 
the unperturbed trajectory must be smaller than 1, in other words that our approximate semiclassical 
estimate of the atomic wavefunction be essentially the same as in the WKB method. 

For illustration, let us consider the situation of a phase grating modulated along a single direction 
in space with a period a, and let rimax be the maximum diffraction order observed in the experiment. 
Condition (8) then reads: 



(9) 
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Fig. 1. — Perturbed (solid line) versus unperturbed (dashed line) atomic trajectories. At time t = U where 
the atoms exit the interaction region with the dii&action grating, the maximum deviation in momentum (resp. in 
position) between the perturbed and tmperturbed atomic trajectories is Apmax (resp. Avmax)- 



which means that the deviation of the atoms from their unperturbed trajectories must be small compared 
to the grating period divided by the maximum diffraction order. Note in particular that condition (9) is 
stronger than the validity condition of the WKB method, which requires that no caustics or focus points 
appear inside the interaction region between atoms and laser, and which reads: 



^rmax < a (10) 

As it happens, the validity condition (8) of the perturbation method embodies the range of applicability 
of the semiclassical method, provided however that the incident atomic de BrogUe wavelength remains 
small compared to the spatial variation scale of the optical potential. Finally, it is possible to transport 
condition (8) from the spatial to the time domain by noticing that Avmax is of the order of Ap^ax/M 
times the typical interaction time r between the atoms and the Ught grating. In this way, one obtains a 
validity condition equivalent to (8) which reads: 

^^T<n (11) 

To conclude, we point out some important features of the above perturbative method. First, the 
method clearly establishes an analogy between atomic diffraction by a weakly modulated optical po- 
tential, and light diffraction by a thin phase grating. Second, it is valid no matter the form of the 
unperturbed atomic trajectories, and thus identically applies to transmission and reflection gratings. 
Third, as will be shown in the following, it generally leads to analytical expressions of the diffrac- 
tion amplitudes because Si only involves the unperturbed atomic trajectories which are often known 
analyticaUy. 
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2.3. Calculation of the atomic diffraction spectrum. — In a real experiment, atomic 
diffraction is observed at a distance from the grating much larger than its typical spatial modulation 
scale (generally of the order of a micron). As a result, the atomic wavefunction (2) evaluated just at 
the exit of the interaction region with the grating does not properly account for the observed diffraction 
pattern. In fact, it is necessary to propagate the atomic wavefunction in the far-field region for obtaining 
the appropriate diffraction spectrum. Here are summarized the main results of Appendix B, where the 
detailed derivation of the diffraction spectrum is presented. 

For illustration purposes, we consider a two-dimensional geometry (x-z-plane) where atoms are 
incident from z = — oo with momentum pi on a diffraction grating with spatial modulation of period a 
in the x direction (the formalism can be generalized straightforwardly to a three-dimensional situation). 
Because of the periodicity of the grating, the atomic diffraction spectrum displays a discrete pattern. 
The diffraction orders are labeled by an integer number n which denotes the momentum transfer nhq 
(with q = 27r/a) from the grating to the atoms in the x direction. It follows from energy conservation 
during the diffraction process that the atomic momentum p(") = associated with the n-th 

diffraction order reads: 



Px = Pi, x + 'iiiq 



(12) 

2 



where pi^x is the component of the incident atomic momentum along the x axis. As shown in Appendix 
B, the diffraction ampUtude a„ associated with the n-th diffraction order can be evaluated from the 
value of the atomic wavefunction on a surface z = Zi located at the exit of the interaction region 
between the atoms and the Ught grating. More precisely, two situations can be distinguished depending 
on the way the atomic wavefunction is evaluated. 

• Derivation from the WKB wavefunction. — In the case where the atomic wavefunction is evaluated 
following the WKB method (Sect. 2.2.1), the diffraction amplitudes are given by (see Appendix B): 



(n) 
Pi 



V(rf)exp(-^p(")-rf ) (13) 



where rf = [xi, zf) belongs to the ideally infinite surface z = zi andpt^zixi) denotes the z-component 
of the atomic momentum associated with the classical perturbed trajectory r{t) satisfying the boundary 
conditions (4). It is important to note that equation (13) does not merely correspond to the Fourier 
transform of the atomic wavefunction after the interaction with the diffraction grating. A difference 
indeed arises due to the supplementary factor in square brackets, which accounts for the angle of inch- 
nation of the trajectories with respect to the normal of the surface z = Zf. 

• Derivation from the perturbation method. — In the case where the atomic wavefunction is evaluated 
following the perturbation method described in section 2.2.2, it is consistent to evaluate the diffraction 
amplitudes by substituting condition (1 1) into equation (13). As shown in Appendix B, this yields: 

a„ = ^^ da;f^(rf)exp (^-^p(") Tf^ (14) 



which is the simple Fourier transform of the atomic wavefunction after the atom-grating interaction. 



N°ll 



ATOMIC DIFFRACTION BY A THIN PHASE GRATING 



1961 



3. Illustration example: diffraction by a Gaussian standing wave. 



We illustrate in this section the capabilities of our perturbation method by considering the well-known 
case of atomic diffraction by a Gaussian standing wave, leading to the nearly-resonant Kapitza-Dirac 
effect. We show that our method recovers the analytical expression for the diffraction orders [1,2] 
at any incident angle, and we prove its equivalence to the Raman-Nath approximation. The validity 
condition (8) of the method is verified by comparison with numerical WKB calculations, and we show 
that the accepted validity conditon of the Ranam-Nath approximation can lead to inaccuracies in the 
evaluation of the phases of the diffraction amplitudes. 



3.1. Calculation of the diffraction spectrum. — Consider a two-dimensional geometry 
where an atomic beam of momentum pi crosses the waist w of a Gaussian standing wave at the inci- 
dence angle 6* (6* = at normal incidence), and assume that the atomic kinetic energy is much larger 
than the height of the optical potential provided by the light field. The Lagrangian describing the atom 
dynamics thus takes the form (1) with: 

Lo = ^Mr2 (15a) 

sV{r) = -^Viexp{-2z'^/w^){l + cos2kx) (15b) 
V 27r 

where k = 27r/A is the wavevector associated with the laser wavelength A. As described in section 
2.2.2, the perturbation method is based on the evaluation of the action integral Si (Eq. (6b)) describing 
the phase-shift undergone by the atoms along their unperturbed trajectories ro(t) = {xo{t), zo{t)) 
which read: 




(16) 

The physical origin of the phase-shift can be visualized in figure 2a where we have projected the tra- 
jectories (16) onto the perturbation potenial (15b). More precisely, is evaluated by substituting 
equations (15b) and (16) into (6b). One thus gets: 

Si = -^ViT{l+(3KD{0)cos2kxO (17) 

where r = Mw/p; ^ is the typical interaction time between the atoms and the laser wave, and Pkd{()) 
is a real parameter describing the amplitude of the phase modulation which depends on the incidence 
angle 6* [11]: 

I3kd{0) = exp (^-^{kw tanef^ (18) 

Note that the Kapitza-Dirac incidence factor Pk d decreases exponentially with the parameter kw tan 9 = 
k{pi^x/M)T, which corresponds to the dimensionless displacement of the atoms along the standing 
wave direction during the typical interaction time r. In particular, /3kd = 1 (maximum value) at nor- 
mal incidence, whereas Pkd tends toward zero at grazing incidence (because of the spatial averaging 
of the potential modulation). 
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Fig. 2. — Unperturbed atomic trajectories (bold lines) projected onto the perturbation potential cl^(mesh). The 
atomic diffraction process results from the spatially modulated phase-shift undergone by the atoms along their un- 
perturbed trajectories where they interact with the perturbation potential, (a) Case of a Gaussian standing wave, (b) 
Case of an evanescent wave reflection grating (x — z coordinates). Note that the potential surface only corresponds 
to the spatially modulated part of the optical potential, and thus takes both positive and negatibe values, (c) Same 
as figure 2b, but a-time coordinates. Note the analogy with the case of figure 2a. 
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Finally, using equations (17), (5) and (2), one obtains: 

tp{xf, Zf) oc exp 



— — eVir/3if£>(^) cos[2fc(xf — 2;f tan^)] 



from which the populations of the diffraction orders follow straightforwardly (see Eq. (14)): 

eViT 



Pkd{0)- 



2h 



(19) 



(20) 
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angle of incidence (mrad) 
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7l/4 
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Fig. 3. — Populations of the diffraction orders n — 1-3 versus the incidence angle (diffraction orders n and 
— n are equally populated), (a) Case of a Gaussian standing wave transmission grating, with w = lOOA and 
sVit/ {2Ji) = 3. The diffraction pattern displays a dramatic dependence on 9 on the scale of a few mrad. (b) By 
contrast, the evanescent wave reflection grating leads to a diffraction spectrum which varies much more smoothly 
with the incidence angle (« = q, £pi,z/ {hn) = 3). 



1964 



JOURNAL DE PHYSIQUE II 



N°ll 



with Jn the n-th Bessel function of integer order. Note that equation (20) exactly corresponds to the 
result derived by Martin et al. [12] using the Raman-Nath approximation. Note also that the estimate 
(20) is subject to the condition (1 1), which can be written in the form: 

^nl^axERT « n (21) 

where Er = h'^k'^/2M is the one-photon recoil energy, and Timax is the maximum diffraction order, 
given approximately by Umax ~ Pkd^^- 

We have represented in figure 3a the dependence of the populations of the three first diffraction orders 
as a function of the angle of incidence 6, for a typical experimental value of the laser waist w = 100 A 
[12]. As reported in [12], the diffraction spectrum is found to display a dramatic sensitivity on the 
incidence angle on the scale of 5 mrad. This property results from the large value of w used in figure 
3a as well as in the experiment of reference [12] (see Eqs. (18) and (20). 

3.2. Comparison with the Raman-Nath approximation. — As shown above, the perturbation 
method allows us to recover exactly the analytical expression for the populations of the diffraction 
orders obtained using the Raman-Nath approximation. We show here that both approaches are actually 
equivalent in the case of a Gaussian standing wave diffraction grating. 

Let us briefly recall the main features of the Raman-Nath approximation [2, 12]. This approach is 
based on the Schrodinger equation describing the interaction between fast atoms of incident momentum 
Pi and the laser standing wave. In the interaction picture and after adiabatic eUmination of the excited 
state, this equation reads: 

indMx,t) = (^-^^+eV[rom^ ^{x,t) (22) 

where V is defined as in equation (15b) and ro(i) as in equation (16). The atomic wavefunction is then 
Fourier-expanded as tpix, t) = a„(t)e^'"'^^, and the diffraction amplitudes a„ are found to satisfy: 



jd„ = — =yiexp(-2tVT^) X (23) 



Finally, the kinetic energy term in equation (23) is neglected compared to the average atom-field cou- 
pUng, provided that the Raman-Nath condition is fulfilled: 

4n^„^£;ij <C e^kdVx (24) 

Hence expression (20) is readily obtained. 

In fact, the Raman-Nath approximation is equivalent to the Schrodinger equation (22) -without the 
kinetic energy term, the solution of which: 

V'(a;f,if) (xexp (^-Je^^ dtV[ro{t)]^ (25) 

is rigorously equivalent to the result (19) obtained using the perturbation method. It is important to 
note however that despite the mathematical equivalence between the Raman-Nath approximation and 
the perturbation method, both approaches are associated with different validity conditions. Indeed, 
whereas condition (21) apphes to the latter, the Raman-Nath approach is associated with equation (24) 
which can be rewritten using the expression of Umax in the form: 
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2nmaxERT «: n (26) 

Thus condition (21) is more severe than condition (26) by a factor of the order of Umax, which can be 
significantly larger than 1 . We show in the following section that as far as the populations of the diffrac- 
tion orders are concerned, conditions (21) and (26) are not easily distinguishable, but that equation (21) 
is actually more accurate than (26) for the estimate of the phases of the diffraction amplitudes. 

3.3. Comparison with WKB calculations. — Here we discuss in more details the validity 
of the perturbation method, by comparing it to the WKB treatment of atomic diffraction. In view of 
section 2, one can express the WKB wavefunction V^vkb in terms of the perturbed wavefunction '0pert 
as: 

^WKB = Vperte'*^ (27) 

where in the Umit of small deviations between the two approaches, 

5(f)ix{) w -r;sin2[2A:(a;f - tan6')]. (28) 

Here rj = ^n^g^^E / % is the small parameter from condition (21). As a result, from equation (27) 
and neglecting possible asymmetries in the diffraction, the diffraction spectrum derived from the pertur- 
bation method appears as the convolution of the WKB diffraction spectrum with the Fourier transform 
of the function: 

exp(i^0) = e"'''/^ ^ i"J„(77/2)exp[4imA;(a;f-^ftan6i)] (29) 

m=— oo 

which only displays even diffraction orders, and has a typical width 5m. k, rj. 

In the hmit rj ^ 1 where condition (21) is not fulfilled, the convolution spectrum is large and thus 
the diffraction spectrum derived using the perturbation method differs significantly from the WKB 
spectrum. In contrast, in the limit <C 1 where the validity condition (21) is fulfilled, the convolution 
spectrum is narrow, and hence the perturbation method is a good approximation. Quantitatively, the 
second order expansion of function (29) reads: 



exp{i6(j)) « (30) 



e 



2 2 
TJ 71 71 

1 -h «— cos[4A;(xf — Z{ tan^)] — — cos[8A;(a;f — tan^)] — — 



As far as the phases of the diffraction amplitudes are concerned, the difference between the WKB and 
the perturbation method thus reduces to a global phase-factor e"*/^. It is also possible to compare the 
populations of the diffraction orders using equation (30). After a straightforward calculation, one finds: 



On.pcrtP — |a„,WKBp ~ ~ Tr + Jn-4 + Jn+4) — {Jn-2 + >/n+2)^] (31) 



t 

16 

where the argument of the Bessel functions is the same as in equation (20). Since the term in square 
brackets is typically of the order of unity, we see that the populations of the diffraction orders derived 
from the WKB and the perturbation method only differ by a fraction of rf . It is important to note that the 
perturbation method, and hence the Raman-Nath approximation, is therefore much more accurate for 
the evaluation of the modulus square than for the estimate of the phases of the diffraction amplitudes. 



1966 



JOURNAL DE PHYSIQUE II 



N°ll 



We have performed numerical WKB calculations of the diffraction spectrum, and compared the 
results with those of the perturbation method. We have confirmed that in the case where condition 
(21) was fulfilled, the phases of the diffraction amplitudes actually differ by the amount —ri/2. The 
populations of the diffraction orders were not clearly distinguishable as expected from the smallness of 
rf. We have been particularly interested in the range of parameters where condition (26) was fulfilled, 
whereas condition (21) was not. In the parameter space that we have explored, we have observed that 
the populations of the diffraction orders were still given to a good approximation by equation (20), but 
that the phases of the diffraction amplitudes actually differed significantly from the WKB estimates. 

In conclusion, it turns out that the vaUdity of both the perturbation method and the Raman-Nath 
approximation are subject to one condition for the phases (Eq. (21)) and to another for the modulus 
square (Eq. (26)) of the diffraction amplitudes. 



4. Diffraction by an evanescent wave reflection grating. 

We consider in this section the atomic diffracton process associated with the reflection grating provided 
by an evanescent wave having a small standing wave component. Similar to the analogy between the 
Kapitza-Dirac effect and Ught diffraction by an acoustic wave, one might expect the evanescent wave 
reflection grating to be similar to the light diffraction grating produced by surface acoustic waves (where 
the periodic undulations of the free surface act as a surface grating). In fact, important differences arise 
between atom and conventional optics concerning the reflection gratings. First, the repulsive potentials 
of atomic mirrors typically vary on the scale of the optical wavelength, which is generally much larger 
than the de Broglie wavelength of the incident atoms. The situation is reversed in conventional optics 
where metalUc or dielectric surfaces achieve spatial changes in the refractive index on a spatial scale 
much smaller than the optical wavelength. Second, from a geometrical optics point of view, atomic 
trajectories display properties different from those of light rays because of the nonzero atomic mass. 
For instance, the fact that an atom can be decelerated until its velocity is zero has no equivalent in 
conventional optics. We show here that in the semiclassical regime of reflection, the evanescent wave 
reflection grating is actually much more closely analogous to the transmission than to the reflection 
grating of conventional optics. More precisely, we prove that it is equivalent to the transmission grating 
produced by a standing laser wave having an Eckart profile (a sech^). Analytical expressions for 
the populations of the diffraction orders at any incidence angle are obtained, and the experimental 
conditions for observing the diffraction spectrum in the thin phase grating limit are briefly discussed. 



4.1. Calculation of the diffraction spectrum. — Consider a two-dimensional geometry 
where an ensemble of laser-cooled atoms of momentum pi is incident on an evanescent wave reflection 
grating [3, 4] having a smaU standing wave component. Let 9 denote the angle between pi and the 
normal of the mirror, and assume that the optical potential height provided by the evanescent wave is 
larger than the incident kinetic energy of the atoms. The Lagrangian describing the atom dynamics thus 
takes the form (1) with [13]: 



eV{T) = eFie-2«^cos2ga; 



(32a) 
(32b) 
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where the wavevectors k and q are typically of the order of the vacuum wavevector fc = 27r/ A associated 
with the laser wavelength A. Following the semiclassical perturbation method described in section 
2.2.2, we first derive the unperturbed atomic trajectories ro{t) = {xo{t) , zo{t)) which are given by 
[14]: 



Xo{t) 

zo{t) 



2k 



Pi 



2MVi 



sech^ (t/r) 



(33) 



where r = M/ npi^z it the typical reflection time of the atom on the evanescent wave mirror. Second, we 

evaluate the action integral 5i (Eq. (6b)) by time-integration of the perturbation potential experienced 
by the atoms along their unperturbed trajectories (see Fig. 2). Using equations (32b) and (33), this 
potential is found to be: 



eVit) = cos[2g(xi + pi,,/Mt)]sech2(t/r). (34) 

By comparison with equations (15b) and (16), it clearly appears that equation (34) is analogous to the 
perturbation potential experienced by the atoms in the Kapitza-Dirac geometry discussed in section 3 
with the foUowing substitutions: 



' 14- 




k- 


' q 


w - 


1/k 



(35) 



exp(— sech (kz). 



This shows that in the limit of a small standing wave component, the evanescent wave reflection grating 
actually behaves as a transmission grating (see Fig. 2). This is because the optical potential associated 
with the evanescent wave varies very smoothly on the scale of the incident de Broglie wavelenght [15]. 
Substituting equation (34) into expression (6b), one readily obtains: 



= — Pew (0)^-^ cos 2qxi 

where Pew{&) is the evanescent wave incidence parameter analogous to equation (18) given by: 

7rtan^g/K 



PEw{e) = - 



sinh(7rtan^g/K) 



(36) 



(37) 



which has a similar asymptotic behaviour as Pkd (Eq.(18)) as 6 tends toward or -k/2. FinaUy, using 
equations (36), (5), (2) and (14), it is straightforward to derive the populations of the diffraction orders 
which are found to read: 



|a„|2 = Jl (ePEwie f-^) . 



(38) 



Note that equation (38) is analogous to expression (20) with the replacements (35), as expected from 
the analogy with the Kapitza-Dirac effect. Note also that the argument of the Bessel functions (or 
equivalently the number of observable diffraction orders) is here proportional to the incident momentum 
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in the ^ direction in the case of the reflection grating, whereas it is inversely proportional to it in the 
case of the transmission grating (see Eq. (20)). This is because in the case of the transmission grating, 
the height of the perturbation potential is fixed and therefore the atomic phase-shift decreases with 
incident momentum because of the 1 /pi^^-dependence of the interaction time r. By contrast, in the 
case of the reflection grating, the height of the effective perturbation potential experienced by the atoms 
is proportional to the square of the incident momentum, as shown by equation (35). For increasing 
momentum the increase of the potential height therefore overcomes the decrease of the interaction 
time, hence the phase-shift increases with pi^. 

Similarly to the case of section 3, the accuracy of equation (38) is subject to condition (21) with 
iT-max ~ s:f3Ew{(^)Pi,z/^K, with the restriction however that the semiclassical description of the re- 
flection process be valid, i.e., pi 2 ^ Hk [14]. In fact, even though condition (21) is required for the 
accuracy of the phases of the diffraction ampUtudes, equation (38) remains valid in a broader range of 
parameters given by a condition similar to equation (26). 

We have represented in figure 3b the dependence of the populations of the three first diffracton orders 
as a function of the angle of incidence 9, for typical experimental parameters. In contrast with the case 
of a Gaussian standing wave (Fig. 3(a)), the diffraction spectrum of the evanescent wave reflection 
grating does not exhibit a dramatic sensitivity on the incidence angle. This characteristic arises from the 
difference in the spatial extensions of the gratings along the z direction (iv = lOOA for the ttansmisson 
grating, tz; = 1/k ss A/27r for the reflection grating). 



4.2. Conditions for an experimental realization. — Here we briefly discuss the experi- 
mental conditions required for the observation of atomic diffraction by an evanescent wave grating in 
the regime of small modulation of the optical potential investigated in the preceding section. Let us 
first examine the implications of the validity condition (21). Using the expression of r and assuming 
that K Ri fc, one obtains as a first constraint: 

Pi,z > 2nL,nfc. (39) 



In order to observe 5 diffraction orders, one thus has to achieve typically pi^z ~ WOfiK. Second, one 
has to take care of the absence of spontaneous emission events during the reflection process. As shown 
in [16], the spontaneous emission probabihty Pgp per reflection reads: 

Psp = ^f^ (40) 

A TiK 

where F is the natural width of the excited state and A is the laser frequency detuning from resonance. 
In order to avoid spontaneous emission, one thus needs a laser detuning A w lO^F. Third, it is 
necessary to realize a sufficiently high optical potential barrier for reflecting the atoms, i.e., Vi w 
lO^Eji. Because Vi is inversely proportional to the detuning A, this requires a large laser intensity, 
typicaUy of the order of 10^ W/mm^ (case of rubidium atoms). This shows that such an experiment 
requires amplification techniques of the evanescent wave using either surface plasmons [17] or thin 
dielectric waveguides [18]. 



5. Application of the method to other experimental situations. 

We have shown so far that the semiclassical perturbation method was a convenient tool for describing 
atomic diffraction by a weakly spatially modulated potential. We mention here possible extensions of 
the method to other situations of experimental interest. 
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• Atomic interferometry. — In atomic interferometry, it is often necessary to derive atomic phase- 
shifts due to small potentials, associated for example with a gravitational field or the rotation of the 
interferometer. In this situation, the expansion (6) can be used to calculate perturbatively the atomic 
phase-shifts using action integrals along the unperturbed atomic trajectories [19]. 

• Time-modulated potential. — Because the central point of the perturbation method consists in eval- 
uating the time-integral of a perturbation potential experienced by an atom along its unperturbed tra- 
jectory, the method applies straightforwardly to the situation where an atom interacts with an optical 
potential having a sufficiently small time-modulated component. For an illustration in the case of a 
time-modulated evanescent wave optical potential, see reference [20]. 

• Multilevel atoms. — Certain experimental situations arise where the multilevel structure of the atoms 
plays an important role. For example, one can be interested in the atomic diffraction process by a 
Gaussian standing wave saturating the transition between a ground-state and a long-Uved excited state. 
One can also consider a situation where the saturation of the optical transition is negligible, but where 
the Zeeman degeneracy of the ground-state is involved (e.g., an atom interacting with a standing wave 
displaying a polarization gradient). In such cases, the atomic state must be described by a spinor the 
components of which are associated with a given atomic internal state. However, in the case where the 
potential responsible for the mixing of the internal states can be considered as a small perturbation, it 
is possible to generaUze our method and to derive the time-evolution of the atomic spinor by integra- 
tion of the evolution operator associated with the perturbation potential along the unperturbed atomic 
trajectories [21]. 

6. Conclusion. 

We have presented a semiclassical perturbation method which allows one to describe in a simple way 
the interaction between an atom and a potential having a small modulated component. This method 
is the analog in atom optics of the treatment of light interaction with thin phase objects in conventinal 
optics. It generally provides clear physical pictures as well as analytical descriptions of the interac- 
tion process. It consists in evaluating the atomic wavefunction by time-integration of the modulated 
potential along the unperturbed classical trajectories of the atoms. A validity condition of the method 
has been given and illustrated on the well-known Kapitza-Dirac effect, where the method has been 
proved to be equivalent to the Raman-Nath approximation. We have also used the perturbation method 
for deriving for the first time an analytical expression for the populations of the diffraction orders of 
an evanescent wave reflection grating at any incidence angle. This perturbation method should prove 
interesting in a broader range of experimental situations, as shown by its application to the treatment of 
atom interaction with time-modulated potentials [20] . 
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Appendix A. 

Semiclassical perturbative derivation of the atomic wavefunction. 

In this appendix we derive equations (2) and (3) of section 2. 1 for the semiclassical atomic wavefunc- 
tion, as well as the perturbative expansion (6) of the integral (3). 
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A.l Semiclassical derivation of the atomic wavefunction. — For a time-independent 
Hamiltonian system, the quantum propagator can be represented by a Feynman path integral [22]. 
Thus, if at time i = , the atom is described by a plane wave of momentum pi: 



ip{r,ti) =exp-(pi -r) 
then at time t = ti,thc atomic wavefunction is given by the convolution product: 



V'(rf , tf) = J driK{rf, if [r;, U) exp ^^p; • r;^ 



(A.1) 



(A.2) 



with the Feynman propagator defined as the path-integral: 



P[r(i)]exp- ( / L{vit),r{t))dt 



(A.3) 



The measure I'[r(i)] signifies that the integration is to be taken over all trajectories r{t) satisfying the 
boundary conditions: 



ri 



(A.4) 



The semiclassical version of the quantum propagator (A.2) arises from a stationary-phase approxima- 
tion. The path integral is then dominated by contributions from classical trajectories since these render 
the phase of the integrand stationary. The phase of the atomic wavefunction (times h) is therefore given 
by the value of the generahzed action: 

S[r{t), ri] = Pi . r; + ' L(r(i), r{t))dt (A.5) 
for the specific initial point ri and the trajectory r{t) which fuUfiU the stationary-phase condition: 



SL 

= SS = Pi ■ 5ri + Sr ■ — 
or 



dr dt dr 



(A.6) 



for any small deviation {6ri,6r{t)) from the path of stationary phase. Because the boundary conditions 
(A.4) impose the relation: 



it follows from (A.6) that: 



6r{ti) = Svi 



P(*i) = ^(ii) = Pi 

dL d aL _ ^ 
, dr dt dr 



(A.7) 



(A.8) 



As expected, one finds that the stationary-phase trajectory r{t) is a classical trajectory (it is solution of 
the Euler-Lagrange equations of motion) which satisfies the boundary conditions (4). The semiclassical 
atomic wavefunction is hence given by equations (2) and (3) of section 2.1. 
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A. 2 Perturbative expansion of the ACTION INTEGRAL. — We consider the method described 
in section 2.2 for deriving the action integral (3) by an expansion in powers of the small parameter e. 
In fact, only the second-order term is non-trivial: the zero-th order term Sq corresponds to the absence 
of the modulated potential and hence to expression (6a), and equation (6b) arises from the fact that the 
term linear in e involving the nonperturbed Lagrangian Lq vanishes by virtue of the stationary phase- 
condition (A. 6) for the unperturbed trajectory ro(<). In order to calculate the second-order term (6c), 
we have to expand both the boundary conditions (4) and the Euler-Lagrange equations of motion for 
the perturbed trajectory up to first order in e. This yields the relations: 



= (r^.^ + n.— ^ ^ 

dr dr J dr 







d \ dLn 



d_ _ 
dr dr J dr 



r=ro, r= 


=rQ, t=to 




dV 




d 


dr 




dr 



(A.9) 



d_\ dLo 
dr J dr 



r=ro, r=ro 



The second-order term of the action is equal to: 

^2=Pi- 



, , /■*'■,/ dLo . dLo dV 1 / d . d \W ^ ^ 



10) 



Note that all the quantities inside the integral are evaluated along the unperturbed trajectory. Finally, 
using (A.9), the Euler-Lagrange equation for ro{t) and integrations by parts, equation (A.IO) is readily 
simphfied to get the result (6c). 



Appendix B. 

Calculation of the diffraction spectrum. 

In this appendix we derive the expressions of the atomic diffraction amplitudes (13) and (14). 

B.l Quantum-mechanical integral theorem of Helmholtz and Kirchhoff. — In a 
potential-free region of space, the wavefunction of an atom of kinetic energy E satisfies the Schrodinger 
equation: 

, 2ME 

VV + — ^V' = 0. (B.l) 

Because equation (B.l) has the same form as the Helmholtz equation in electromagnetic theory, it 
is possible to use a quantum-mechanical version of the integral theorem of Helmholtz and Kirchhoff 
[23] to express the atomic wavefunction in the far-field region (which describes the diffraction pattern) 
from its value on a boundary surface E located in the free-field region reached by the atoms after their 
interaction with the diffractron grating. One thus has: 

V'(r)farfieid = ^JdEn- [GE{r, rf )VV'(rf ) - V(rf ) VG^Cr, rf )] (B.2) 

where n denotes the outward-pointing normal to the surface S, and rf corresponds to the endpoints of 

the classical trajectories along which the action integral (3) is calculated. The integral in equation (B.2) 
involves on one hand the quantum propagator (r , r f ) for a free particle of kinetic energy E, solution 
of: 
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2M _ 



-E 



Gsir, rf ) = -47r(5(r - rf ) 



(B.3) 



(S is the Dirac delta function), and on the other hand the value of the atomic wavefunction and its 
spatial derivative n • Vip on the surface S. More precisely, is evaluated either from the WKB or 
the perturbation method using the relation: 



Vij{rf) = -pf(rf)i/'(rf) (B.4) 

where Pf (rf ) is the local atomic momentum at the endpoint point rf of the classical trajectory. Note 
that whereas the direction of Pf (rf ) generally depends on rf, its modulus is constant, equal to that of 
the incident atomic momentum pi (because of energy conservation in the diffraction process). 



B.2 Diffraction by a two-dimensional periodic grating. — We consider here the partic- 
ular case of a two-dimensional diffraction grating located in the xOz plane (or equivalently a three- 
dimensional grating with translational invariance along the y-axis). We assume that the grating has a 
finite extent in the z-direction, but is infinite in the .x-direction along which it is spatially modulated 
with the periodicity a. Using the notations of section 2, one thus has: 

Va;, V{x + a) = V{x) (B.5) 
It follows from equation (B.5) and the Bloch theorem that the atomic wavefunction will take the form: 

Vx, z ip{x + a,z) = e'P''«"/^^/>(x, z) (B.6) 

with pi^x the component of the incident atomic momentum along the x-axis [24]. In order to take 
advantage if this property in the derivation of the diffraction spectrum, it is convenient to define the 

surface E of equation (B.2) as a line z = zs (n = —e^ is thus independent of rf) and to express the 
two-dimensional quantum propagator Ge{y, rf) in the free-field region {z > Zf) in the form [23]: 

G£;(r,rf)=i J dp^,^ exp (^^p ■ {r - ri)j (B.7) 

where Pz = 2ME — p^ (the integration is restricted to the interval \px \ < \/1ME because the 
evanescent wave components of the proagator do not contribute to the far-field wavefunction). Using 
equations (B.4), (B.6) and (B.7), and dividing the integration range of equation (B.2) into intervals of 
length a, one obtains: 



V'(r)farfield = V - / dxf / Ap, 



Pz 



ilj{X{, Z{) X 



X exp ( ^n(pi,x - Px)a] exp f ^p • (r - rf ) ) . (B.8) 
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Using the relation: 



+00 



+00 



^ exp Tn{Pi,x - Px)a = %q ^ 5{jPx - Pi,x - nhq) 



(B.9) 



n— — 00 



and expression (12) for the momenta associated with the different diffraction orders, equation (B.8) 
finally yields: 



+00 



V'(r)farfield = ^ On exp ( tCP^"-" ' 1") 



(B.IO) 



where o„ is the diffraction amplitude associated with the n-th diffraction order, which reads: 



1 r 



1 + 



Pz 



^/'(rf ) exp 



_lp(n) . 

h 



rf 



(B.ll) 



Even though equation (B.ll) may be used to derive the atomic diffraction spectrum in the framework 
of the perturbation method, it seems more consistent to retain only the terms of (B.ll) which correspond 
to the accuracy range of the method. As previously discussed in section 2.2.2 (Eq.(l 1)), in the validity 
domain of our semiclassical method one has: 



^fL^ « a « 1 



and hence: 



which yields: 



Pi 



< 1 



Pf,z 

(n) 
Pz 



Et 



Ap: 



,2 

max 



Pt 



< 1 



Ap2 

-t max 



<c 1 



(B.12) 



(B.13) 



(B.14) 



Finally, one obtains the expression of the diffraction amplitude for the atomic wavefunction derived 
using the perturbation method: 



an = ^J da;f^(rf) exp (^"^P^"^ • ^f^ 



(B.15) 



which is merely the Fourier transform of the wavefunction evaluated on the line z = Zi, after interaction 
with the diffraction grating. 
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